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ABSTRACT 

We analyze a high resolution spectropolarimetric dataset collected for the He-weak 
B3p IV star HR 2949. The Zeeman effect is visible in the circularly polarized com¬ 
ponent of numerous spectral lines. The longitudinal magnetic field varies between 
approximately —650 and -1-150 G. The polar strength of the surface magnetic dipole is 
calculated to be 2.4ta2 ^G. The star has strong overabundances of Fe-peak elements, 
along with extremely strong overabundances of rare-earth elements; however, He, Al, 
and S are underabundant. This implies that HR 2949 is a chemically peculiar star. 
Variability is seen in all photospheric lines, likely due to abundance patches as seen 
in many Ap/Bp stars. Longitudinal magnetic field variations measured from different 
spectral lines yield different results, likely a consequence of uneven sampling of the 
photospheric magnetic field by the abundance patches. Analysis of photometric and 
spectroscopic data for both HR 2949 and its companion star, HR 2948, suggests a 
revision of HR 2949’s fundamental parameters: in particular, it is somewhat larger, 
hotter, and more luminous than previously believed. There is no evidence of optical or 
ultraviolet emission originating in HR 2949’s magnetosphere, despite its moderately 
strong magnetic field and relatively rapid rotation; however, when calculated using 
theoretical and empirical boundaries on the initial rotational velocity, the spindown 
age is compatible with the stellar age. With the extensive phase coverage presented 
here, HR 2949 will make an excellent subject for Zeeman Doppler Imaging. 

Key words: Stars : individual : HR 2949 - Stars: magnetic fields - Stars: binaries: 
visual - Stars: chemically peculiar - Stars: abundances - Stars: early-type. 


1 INTRODUCTION 

Magnetic, massive stars have recently emerged as a ma¬ 
jor new sub-class of early-type stars. The Magnetism in 
Massive Stars (MiMeS) Survey Component has established 
that approximately 10% of OBA -type stars ho s t large - 
scale, organized magnetic fields JCrunhut et al.l l2012bl) . 
These fields range in strength from hundreds of G to sev¬ 
eral kG at the magnetic poles, and are capable of confin¬ 
ing their radiative, ionized winds within stellar magneto- 
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spheres lIud-Doula fc Owockil |2002| ). This both reduces a 
star’s mass-loss rate, and spins it down due to angular 
momentum loss via the co rotating magnetosphere jMestell 
1 19681 : lud-Doula et al.l l2009l l , a process that has been mea- 
sured in the case of the He-stron g B2p stars HD 3777 6 
llMikulasek et al.]l2008h and ct Ori E (iTownsend et al.lboid ) . 
This rapid spindown has consequences for stellar evolu¬ 
tion, in particular strongly inhibiting mixing llMevnet et al.l 

Im3). 

Stars detected as magnetic by the MiMeS Survey Com¬ 
ponent were moved into the Targeted Component, which 
aimed to collect large datasets of high-resolution spectropo- 
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larimetry for all stars of interest. One such star is HR 2949. 
For many years this mid-B type star was considered to be a 
standard star, and no variability was suspect ed in either its 
integr ated light or spectral features. However. I Rivinius et al.l 
(1200311 pointed out that the star is actually variable both 
photometrically and spectroscopically, suggesting it to be a 
He-variable and/or a He-weak star. 

The He-variable stars, b oth He-weak a nd He-strong, are 
almost invariably magnetic l|Prestonlll974h . and indeed the 
dis covery of a magnetic field in HR 2949 was announced 
by iRivinius et al.l ( 201ir) . who noted the apparently non- 
sinusoidal variation of the longitudinal magnetic field curve, 
along with correlated He line profile variations, implying 
chemical abundance patches as commonly seen in the mag¬ 
netic Ap/Bp stars. 

Many magnetic, massive stars display variable line emis¬ 
sion in optical and/or ultraviolet resonance lines . This is an 
obser vational signature of their magnetospheres (IPetit et al.l 
I 2 OI 3 II . Such emission is not universal amongst the magnetic 
stars, and HR 2949 is amongst the stars that do not dis¬ 
play it. Determining the conditions under which magnetic 
stars do and do not show emission requires precise determi¬ 
nation of the stellar, rotational, and magnetic properties of 
the individual stars. 

In this paper we analyze the photometric, spectropolari- 
metric, and spectral data available for HR 2949 in detail. As 
a control we perform a simultaneous analysis of HR 2948, 
a visual binary companion of similar spectral type but with 
no detectable magnetic field. The stars are of similar magni¬ 
tudes {V = 4.53 and 4.78, respectively), and are separated 
by 7.3”. 

In Section 2, we describe the observations. In Section 

3, we define a rotational ephemeris for HR 2949. In Section 

4, we derive stellar parameters for HR 2949 and its binary 
companion, HR 2948, along with the inclination angle of 
HR 2949’s rotational axis. HR 2949’s spectral variability is 
examined in Section 5. Mean surface abundances for both 
stars are presented in Section 6. The measurement and anal¬ 
ysis of the magnetic data is presented in Section 7, where 
we analyze HR 2949’s magnetic field using both individual 
spectral lines as well as single-element least-squares decon¬ 
volution profiles, and then present a model for the dipolar 
component of the photospheric magnetic field. In the discus¬ 
sion in Section 8 we consider the implications of an appar¬ 
ent discrepancy in the ages of the two stars; interpret HR 
2949’s spectral and photometric variability in the context of 
the star’s magnetic field; and explore the theoretical charac¬ 
teristics of HR 2949’s magnetosphere, including the extent 
of magnetic confinement, degree of rotational support, and 
estimated spindown time. Finally, our conclusions are sum¬ 
marized in Section 9. 


2 OBSERVATIONS 

Spectropolarimetric observations of HR 2948 and HR 
2949 have been obtained with ESPaDOnS at the 3.6m 
Canada-France-Hawaii Telescope (CFHT) in the 2009/2010, 
2012/2013, and the 2013/2014 fall/winter seasons. For each 
observation, the instrument delivered four Stokes I (in¬ 
tensity) spectra and one Stokes V (circular polarization) 
spectrum from 370 to 1040 nm, with a resolving power of 


Table 1. Log of ESPaDOnS observations. When binned spectra 
have been used in the analysis, HJD and phase are given as a 
range from the first to the last spectrum. Exposure times corre¬ 
spond to the total exposure time of all polarized sub-exposures. 
Peak signal-to-noise ratios (SNRs) are per spectral pixel in the 
binned spectra. 


Obs. 

Date 

HJD 

- 2450000 

Phase 

^exp 

(s) 

Peak 

SNR 


HR 2949 



2010-01-01 

5198.0726 

0.020 

600 

826 

2010-01-03 

5200.0880 

0.076 

600 

1245 

2010-01-25 

5222.0194 

0.566 

600 

1397 

2010-01-26 

5222.8473 

1.000 

600 

1372 

2010-01-26 

5223.0030 

0.082 

600 

1279 

2010-01-27 

5223.9870 

0.597 

600 

1306 

2010-01-28 

5224.9357 

0.094 

600 

1284 

2010-01-29 

5225.9513 

0.626 

600 

1387 

2010-01-30 

5227.0052 

0.178 

600 

1103 

2010-01-31 

5227.9687 

0.683 

600 

1040 

2010-02-01 

5228.8705 

0.156 

600 

1121 

2010-02-02 

5230.0433 

0.770 

600 

1032 

2010-02-26 

5253.9359 

0.288 

600 

1375 

2010-02-27 

5254.9081 

0.797 

600 

1225 

2010-02-28 

5255.9602 

0.348 

600 

1281 

2010-03-03 

5258.9242 

0.901 

600 

758 

2012-09-28 

6199.1307 

0.489 

600 

1510 

2012-11-24 

6256.1651 

0.370 

600 

1298 

2012-11-26 

6258.1399 

- 0.1670 

0.404 

- 0.419 

2400 

1105 

2012-11-30 

6262.1494 

- 0.1683 

0.505 

- 0.515 

1800 

794 

2014-02-15 

6703.9418 

- 0.9746 

0.966 

- 0.983 

2240 

1000 

2014-02-16 

6704.8340 

- 0.8779 

0.434 

- 0.457 

2880 

2256 


HR 2948 




2010-10-19 5489.1147 

_ 

1800 

1540 

- 0.1329 

- 




Table 2. Log of FEROS observations of HR 2949. Spectra dis¬ 
carded from the analysis are marked with an asterisk next to the 
calendar date. 


Obs. 

Date 

HJD 

-2450000 

Phase 

Exp. 
Time (s) 

1998-04-12 

1151.7704 

0.150 

240 

1999-01-18 

1196.7644 

0.723 

300 

2003-02-16* 

2686.5113 

0.206 

300 

2003-02-17* 

2687.5075 

0.728 

300 

2004-05-03* 

3129.4594 

0.268 

360 

2004-05-03 

3129.4646 

0.270 

240 

2004-05-03 

3129.4680 

0.272 

240 

2006-02-08 

3774.6919 

0.307 

180 

2006-02-08 

3774.6946 

0.308 

180 
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AA/A ~ 65, 000 at 500 nm. Most of the observations (25) of 
HR 2949 were acquired in the context of the Magnetism in 
Massive Stars (MiMeS) Large Program; the remainder (15) 
were acquired in the context of the BRITE-Constellation 
Spectropolarimetric Snrvey (Pis C. Neiner and G. A. Wade; 
Neiner et ah, in prep.) In total 40 observations of HR 2949 
and 3 of HR 2948 are available. In order to increase the 
SNR, observations acquired on the same night were binned, 
yielding 22 measurements of HR 2949 and 1 of HR 2948 
(the maximum time spanned by a given sequence of binned 
spectra is approximately 1 hour). The binned spectra are 
summarized in the log of ESPaDOnS observations in Ta¬ 
ble [T] 

A total of nine spectra of HR 2949 have been obtained 
with the EEROS spectrograph at La Silla providing a resolv¬ 
ing power of -^48 000 over a spectral range of 375 to 890 nm. 
Four spectra, obtained with EEROS mounted at the ESO 
1.52m and MPG 2.2m telescopes, were already described by 
iRivinius et al.l ll2003f l. Subsequent spectra were obtained at 
three more epochs with EEROS mounted at the MPG 2.2 
m telescope. These data are available from the ESO public 
archive. Three spectra (those acquired in 2003, and the first 
acquired in 2004) were overexposed and excluded from the 
analysis. The log of EEROS observations is given in Table[21 


3 EPHEMERIS 

To constrain HR 2949’s rotational period two independent 
data sets were used. First, there is the Hipparcos pho- 
tometry, for which a period of 1.9093 d was published by 
iKoen k. Ever! ll2002tl . The Hipparcos photometry is shown 
in Fig. [U phased with the spectroscopic ephemeris found 
below. The mean Hp magnitude, 3.75 mag, is much brighter 
than the individual V magnitudes of the two stars (4.53 and 
4.78 for HR 2949 and 2948, respectively), and is consistent 
with both stars having been observed by Hipparcos simul¬ 
taneously. Our ESPaDOnS observations of HR 2948 span 
only 26 min, so variability of this star cannoted be ruled 
out. However, as explored in detail in the remainder of this 
paper, the spectroscopic and magnetic variations of HR 2949 
follow the same period as the photometric variation. There¬ 
fore we suspect that the photometric variability is entirely 
due to HR 2949. 

The data in Fig. [T]have been cleaned by first discarding 
all measurements with Hp > 3.77, as in such cases one of the 
two stars was not fully in the aperature. Re-analysis of the 
cle aned data, ba s ed up on Scargle’s statistics as implemented 
by iKaufer et al.l lll99d L gives P = 1.9083 ± 0.00025 d. For 
presentation purposes, the data were then sigma-clipped by 
removing observations with uncertainties greater than twice 
the mean, and further cleaned by removing points lying more 
than 2a from a 2"‘^-order sinusoidal fit to the remaining 
points. 

In the spectra themselves, the He lines show the 
strongest variability. The ESPaDOnS and EEROS spectra 
were combined into a single set, and the equivalent width 
(EW) of several lines measured, namely He l 400.9 nm and 
402.6 nm together as combined EW, as well as He l 438.8 nm 
and 447.1 nm individually. As the spectra are relatively 
sparse across a long time baseline (1999 - 2009), severe 
aliasing makes it impossible to pick out a single peak in 


3.72 

3.73 
o 3.74 

3.75 

3.76 

3.77 


-0.5 0.0 0.5 1.0 1.5 

Phase (P = 1.90871 d) 

Figure 1. Hipparcos photometry (black circles) phased with the 
ephemeris in equation[T] The solid red line indicates the 2"'‘^-order 
sinusoid used to tr-clip the data. 

the periodogram based on these data alone. However, the 
aliases are relatively narrow, and there is only one peak, at 
P — 1.90871 ± 0.00007 d, that is in reasonable agreement 
with the Hipparcos period. We adopt this period, being the 
most precise. 

As epoch, we choose the ESPaDOnS data showing the 
strongest positive longitudinal magnetic field, so that the 
ephemeris is 

Emax(B)(HJD) = 2455222.84726 -P 1.90871(7) x E. (1) 


4 STELLAR PARAMETERS 

HR 2 949 is usuall y classihed as a B5 or B6 dwarf to subgiant 
star llHouk|[l982l ). However, magnetic, chemically peculiar 
Ap/Bp stars are not well-suited to the usual calibrations for 
spectral type, as the principle of spectral classihcation (to 
strictly classify a star by nearest resemblance of the spec¬ 
trum to a standard star) yields misleading results when con¬ 
fronted with highly over- or under-abundant atomic species, 
especially when lines commonly used for spectral classihca¬ 
tion (such as He) are affected by chemical peculiarity. This 
necessitates a more careful analysis to constrain the physical 
parameters. 

We use the TLUSTY NLTE grid of mo del atmospheres 
for B-stars, BSTAR2006 llLanz fc Hubenvll2007l ). generat¬ 
ing our own grids via linear interpolation between models 
of adjacent Tgff or logg, and convolving the model spec¬ 
tra with the ESPaDOnS instrumental broadening and the 
stellar rotational broadening. As the BSTAR2006 grid does 
not extend to effective temperatures cooler than 15 kK, we 
supplement t his analysis with the LTE grid calculated by 
ICoelhol ll2014h . 

4.1 Effective temperature 

One standard approach to determine Tee is by comparing 
spectral energy distributions (SEDs) to photometric data. 
While most published hux measurements give only the com¬ 
bined magnitude and hux of the HR 2948-1-2949 system, 
there is one published uvby dataset cont aining individua l 
measurements for each stellar component llLindroosI Il9^ ) . 
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Table 3. Stellar parameters. The second column gives the ref¬ 
erence. Numerical references point to the section in this work in 
which the parameter is derived. Alphabetical reference s point to 
the literature according to the following reference key: af iLindroos 


1983), b) iNieval l(2o3)i c) Ivan LeeuwenI (l2007t ). d) iHog et ~ 


| 2000 |). 



Ref. 

HR 2949 

HR 2948 

V (mag) 

a 

4.53 

4.78 

E(B - V) 

4.1 

0.04 

BC (mag) 

b 

-1.80±0.06 

-1.05±0.21 

TT (mas) 

c 

7.18±1.06 

d (pc) 

c 

139 

00 

+ 1 

PM (RA) (mas yr~^) 

d 

-19.4±1.2 

-10.6±1.4 

PM (Dec) (mas yr“^) 

d 

24.1±1.2 

16.3±1.4 

Mv (mag) 

4.2 

^•-*^^- 0.34 

' J - y ^_ 0.34 

Mboi (mag) 

4.2 

- 0.41 

-1 

-^•^^- 0.56 

Teg (kK) 

4.1 

18.5±0.5 

13.6±1.2 

log (Lt/LQ) 

4.2 

3.08/“;i« 

^•oy-0.18 

log 9 (cgs) 

4.2 

4.10 ± 0.15 

4.24±0.15 

Age (Myr) 

4.2 

11 + 15 
-^^-6 

112 / 2 ° 

RARq) 

4.2 

0 q+0-7 

0 'v+1.2 
'^•'-0.6 

M* (Mq) 

4.2 

6.0±0.3 

4.3±0.3 

V sin 2 (kms“^) 

4.3 

61 ± 5 

114 ± 5 

Veq (kmS“l) 

4.3 

90 

-30 

- 

Magnetic Parameters 

i n 

4.3 

40+23 


p n 

7.2 

50±16 

- 

Bp (kG) 

7.2 

2.4±0.3 

- 


These data are reproduced in Fig. [J] Our analysis suggests 
these data were mislabelled in the original catalogue. The 
reasoning is explained below. 

The photometr y was de-reddened using a standard 
ICardelli et al.l |l98^ reddening curve, assuming a small red¬ 
dening of E (B — V) = 0.04 given the stars’ relative proximity 
to the Sun l|van LeeuwenI E oOtI ). 

Chemically peculiar magnetic stars are UV-deficient 
compared to non-magnetic stars of similar spectral type. 
Since the measured values are to be compared to synthetic 
spectra calculated without taking the magnetic field and 
chemical pe culiarities into account, correction factors must 
be applied l|Kochukhov et aLll200,T >. However, as the mag¬ 
netic field strength is moderate (see Section 7), this correc¬ 
tion is minimal for HR 2949 (smaller than the data points 
in Fig. O). 

Using the designations given by iLindroo^ . comparison 
of the photometry to model SEDs suggests a temperature 
of 13.5±0.5 kK for HR 2949, and 18.5±0.5 kK for HR 2948 
The m odel spectral libra ries provided by iLanz &: Hubenvl 
1 I 2 OO 7 I) and lCoelhol ll2014h give essentially identical results. 

An alternate means of determining Teg is via ioniza¬ 
tion ratios. We performed this analysis using ESPaDOnS 
and FEROS observations for a selection of ions (Si ll/Si ill, 
Fe ll/Fe lll, and P ll/P lll), for each case taking the weighted 
mean of measurements from all observations. Equivalent 
Widths (EWs) of the same lines were then measured in 


the synthetic BSTAR2006 spectra, with uncertainties in 
the model EW ratios coming from the range in surface 
gravity (determined below). The observed EW ratios are 
compared to the model EW ratios in Fig. [S] While the 
results differ somewhat for different lines, they indicate 
Tss= 18.5±1.5 kK for HR 2949, clearly inconsistent with 
the photometric results. 

In the course of the abundance analysis (see Section 
6), a similar ionization balance analysis was performed for 
both stars using ATLAS9 and zeeman. This analysis indi¬ 
cated Teff^ 18.3 ± 0.8 kK for HR 2949, and T^s= 13.6 ± 1.2 
kK for HR 2948. These results are, once again, inconsistent 
with the photometric results. 

There may be ionizat i on im balances in the Si H/HI 
lines. iBailev fc Landstreed (l2013l) found abundances deter¬ 
mined using Si H/HI lines to be discordant by up to 1.5 dex 
in all normal B and Bp stars, and concluded that vertical 
stratification is the likely reason. Whether this is the case 
for other elements in late type Bp stars does not seem to 
have been investigated. Furthermore, there are discrepan¬ 
cies in Teff determined spectroscopically and spectrophoto- 
metrically in some other cases. However, the effect is rather 
small : (e.g., 1000 K in the case of HR 7355; iRivinius et al.l 
I2OI3I) and in some cases entirely absent (e.g. HR 5907, for 
which spec tral modeling and spe ctrophotometry give identi¬ 
cal results: fCrunhut et al.|[2012al ). While we did not conduct 
a systematic search for signs of stratification, in the course 
of our abundance analysis, we did not detect any strong 
discrepancies between fits for different ions of the same ele¬ 
ments, suggesting that stratification, if present, is not partic¬ 
ularly pronounced. Furthermore, our photometric and spec¬ 
troscopic analyses yield essentially identical results. Given 
these considerations, we believe that while the assumption 
of ionization balance may well prove to be incorrect, any 
errors thereby introduced into our determination of TeS are 
likely to be on the order of our uncertainties. 

The discrepancy between the photometric and spectro¬ 
scopic determinations of can be easily reconciled if the 
photometric data was mis-labeled. There are reasons to sus¬ 
pect this may be the case: first, the depths of the H Balmer 
lines are clearly inconsistent with the photometric Teff, but 
consistent with the ionization balance Teg; second, bot h 
stars are listed in the 2MASS catalogue (ICutri et alj 200^ 
with HR 2949 the brighter of the two, wh ereas luindroosl 
lll983l) lists it as the dimmer star; finally, iLindroo^ does 
not provide coordinates for his targets, making it somewhat 
more likely for a typographical error to have entered in. 

We thus take HR 2949 to be the hotter, and brighter, 
of the two stars. In light of this, we reclassify HR 2949 as a 
B3p IV star, and HR 2948 as B6 IV. 


4.2 Luminosity, radius, mass, and surface gravity 

To obtain the luminosity, we first calculate absolute magni¬ 
tudes Mv from the apparent V magnitudes and the distance. 
From Mv we calculate the bolometric magnitude Mboi, us- 
ing the bolo metric correction BC from the calibration of 
iNieval (l2013f) . Luminosities are then calculated from Mboi. 
For HR 2949, we find log L/Lq = 3.08toi3; for HR 2948, 
logL/Lo = 2.69t°;?^ 

HR 2949’s position on the Hertzsprung-Russell Dia¬ 
gram (HRD) (Fig. |4]) in comparison with the Geneva evo- 
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HR 2948 



HR 2949 



Figure 2. uvby photometry for HR 2948 (top) and HR 2949 
(bottom) compared to the best-fit model (black), corresponding 
to the parameters in Table [3| As discussed in Section 4.1, the 
original photometric data seems to have been mis-labeled. 



Figure 3. Ionization balances as functions of Tefj for HR 2949. 
Curved lines indicate model EW ratios for \ogg as given in Ta¬ 
ble [Sj dotted lines indicate the EW ratios at the Icr uncertainties 
in \ogg. Filled circles indicate the weighted mean of the mea¬ 
sured EW ratios, open circles the Icr uncertainties in the weighted 
means. The solid vertical line indicates the mean Teff, while the 
vertical dashed lines indicate the standard deviation from the dif¬ 
ferent measurements. 

lutionary tracks and isochrones (jEkstrom et al.|[2012l) im¬ 
plies a zero-age main-sequence mass of 6.0 ± 0.3 M©, a ra¬ 
dius of S.O^Q'g i? 0 , a surface gravity of logy = 4.25lon7 5 
and an age of Myr. HR 2948 should have Mzams = 

4.3±0.3 M©, = 3.7lJ;g i?©, logy= 3 . 941 q; 22 , and an age 

of Myr. This large discrepancy in the apparent ages 

of the two stars is considered further in Section 8 . 

Both stars are in the instability strip for slowly pulsat- 



4.40 4.35 4.30 4.25 4.20 4.15 4.10 

log(^eff) (K) 


Figure 4. Positions of HR 2949 (red circle) and HR 2948 (blue 
square) on the Hertzsprung-Russel diagram. Solid lines indicate 
the evolutionary tracks for 4, 5, and 6 Mq models. Dashed and 
dotted lines indicate isochrones as per the legend. The grey line 
indicates a pre-main sequence isochrone. 

ing B stars. Only HR 2949 can be evaluated for line profile 
variability. We detect no sign of pulsational variability in 
this star: as detailed in section 5, all line profile variability 
can be explained as a consequence of chemical abundance 
spots. 

HR 2949’s surface gravity is consistent with model fits 
to the H/3 line (Fig. [5] bottom), which yield the best fit for 
log^= 4.10 ± 0.05. Care is required in fitting models to H/3 
in ESPaDOnS spectra, as the red wing of H(3 crosses two 
spectral orders. In determining logg, we first merged the 
un-normalized H/3 lines across the two spectra, then normal¬ 
ized the line via a linear fit to nearby continuum regions: this 
avoided warping of the line profile during the order-by-order 
polynomial normalization procedure typically used for ES¬ 
PaDOnS data. As the wings of H/3 are weakly variable (see 
Section 5), we selected a phase found in the middle of the 
range of variation. In the course of this analysis, it was clear 
that the choice of continuum regions for renormalization of 
the line had an effect on logg larger than the formal er¬ 
ror bars. To account for this source of systematic error, we 
increase the uncertainty by 0 . 1 , such that our final surface 
gravity for HR 2949 is log (/= 4.10 ± 0.15. 

The same analysis for HR 2948 yields a higher surface 
gravity, log (/= 4.24 ± 0.15, than is obtained from the HRD. 
However, the results are again compatible within error bars. 

There is some small blue-red asymmetry in the fits to 
the H/3 shown in Fig. [S] In the case of HR 2949 this may 
be a consequence of the weak variability in the line wings, 
however, as the asymmetry affects both stars it is more likely 
that this is an artifact introduced by normalization. 

4.3 Rotation 

The high-resolution ESPaDOnS and FEROS data enable 
very good constraints on usini via line profile modeling. 
This was performed for several species, including Mg ll, 
C II, Si II, and Si lll. Mean line profiles created from all 
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HR 2948 H/S 



HR 2949 HjS 



Figure 5. H/3 lines of HR 2948 (top) and HR 2949 (bottom) 
as compared to model fits. Residuals are shown below the 1- 
dimensional spectra. Vertical black lines demarcate the region 
excluded from the fit. The best-fit models are in red; blue indi¬ 
cate the ±0.15 dex uncertainty. 


HR 2948 Si II 637.1 nm 




Figure 6. Line profile models for the two stars. 


observations were used in order to average out the distor¬ 
tions caused by abundance spots. For HR 2949, we find 
V sini = 58±4 kms“^; for HR 2948, v sini = 112±4 kms“^. 
An independant determination of usini was performed in 
the course of the abundance analysis, obtaining slightly 
higher velocities, although identical within error bars: we 
thus expand our uncertainties to encompass systematic er¬ 
ror, settling on vsini = 61 ± 5 hms”*^ for HR 2949 and 
wsini = 114±5 kms“^ for HR 2948. Fig.[6]shows model line 
profiles spanning these parameters, compared to the mean 
Si II 637.1 nm line. 

In the case of HR 2949, there is some slight disagree¬ 
ment in the wings between the observations and the models. 
This is likely a consequence of saturated local line profiles 
associated with chemical spots in which abundances are in 
excess of the photospheric mean. As discussed below, this 
phenomenon affects essentially all metallic lines in the star’s 
spectrum. 

In the course of obtaining wsini, the systemic velocity 
was determined to be Vsys = 22 ± 2 kms“^, and is identical 
for both stars. 

In the case of HR 2949, as both Prot and usini are 
known, the stellar radius P* is constrained by i. With 
Prot = 1.90871 d, this yields P* sini = 2.18 Pq. The radius 
found above via the HRD would then indicate an inclination 
i = 421^2° I Sind consequently an equatorial rotational ve¬ 
locity Ueq = 91 I 3 Q kms“^. 


5 SPECTRAL VARIABILITY 

HR 2949’s spectrum displays a rich variety of line profile 
variability (LPV), with distinct patterns in different groups 
of chemical elements. Fig. [7] shows dynamic spectra and 
equivalent width (EW) measurements for He, O, Si, and 
H. Dynamic spectra and EWs for C, N, Ne, Mg, P, S, Cl, 
Ca, Ti, Cr, and Fe are shown in the Appendix in Figs. lAll 
and IA2I 

The dynamic spectra are calculated as the residual 
flux of each line as compared to a mean line profile. Ex¬ 
cesses in absorption at a given Doppler velocity appear as 
flux deficits, while deficits of absorption appear as pseudo¬ 
emission. This residual ID flux is plotted in the bottom pan¬ 
els of Fig. [T] The residual flux is then mapped to a grayscale 
colour bar, shown next to the ID flux, with dotted red hor¬ 
izontal lines indicating the mapping from the ID flux to the 
colour bar, such that stronger absorption is darker while 
pseudo-emission is lighter. The color-mapped flux is then 
plotted as a function of rotational phase, with the spectra 
binned into 16 phase-bins, a number chosen so as to reach 
a compromise between higher SNR and phase resolution. 
Black vertical lines represent Tusini. 

The strongest, simplest variations are seen in He (first 
row of Fig. [7|). An absorption excess crosses the stellar disk 
between phases 0.5 and 1.0. Comparison to magnetic data 
(see Section 7) indicates this is associated with the negative 
magnetic pole. O appears to be distributed in a similar fash- 
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Figure 7. Dynamic spectra (left panels) and EWs (right panels) 
of He I, O I, Si II and ill, and H/?. The dynamic spectra are rel¬ 
ative to mean line profiles. The colour-mapping scheme is shown 
relative to the ID residual flux, with the limits of the colour map 
shown with dotted red horizontal lines. The dynamic spectra are 
plotted in order of rotational phase, with the residual flux binned 
into 16 phase bins. Black vertical lines indicate diusinz. 
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Figure 8. Temporal variance spectra for (bottom-top) Ha, H/3, 
and H7. Solid black vertical lines indicate diusinz. Dashed black 
vertical lines indicate the Kepler corotation radius di/?K de¬ 
fined in Section 8.3.1). 


ion (Fig. 0 second row), as do C, Ne, and S (see Fig. lAll in 
the Appendix). A second, weaker absorption excess may be 
present around phases 0.2 to 0.5. 

Si II and Si iii LPV anticorrelates with He (Fig. [7] third 
and fourth rows). Comparison of the EWs of He and Si lines 
makes the anticorrelation quite obvious: peaks and troughs 
of the EW variations are separated by ~0.5 of a rotational 
cycle. Similarly, variation in the core of the H^d line also 
seems to approximately anticorrelate with He (Fig.[71 second 
panel from bottom). 

The most complex variations are seen in N, Mg, Ti, Ca, 
Cr, and Fe. These are shown in the Appendix in Fig. IA2I 
These elements neither precisely correlate nor anti-correlate 
with He. Furthermore, their EWs do not exhibit smooth 
single-wave variations as in He or Si lines. 

The wings of H/3 (Fig. 0 bottom panel) vary in a dif¬ 
ferent way, with red and blue wings moving towards and 
away from the continuum approximately simultaneously. 
This variation is extremely weak. To evaluate the statistical 
significance of the Balmer lin e variability, we creat ed Tem¬ 
poral Variance Spectra (TVS: [Fullerton et al.|[l99^ for Ha, 
H/d, and H 7 (see Fig. A TVS evaluates the deviations in 
individual pixels as compared to scatter in the continuum. 
In order to minimize spurious variability caused by polyno¬ 
mial normalization, in the case of H/I the same procedure 
was followed as described in Section 4 for measnrement of 
logs- 

In all cases, the continnum varies at approximately the 
0.2% level. The strongest variability is seen in the line core. 
This is similar in character to TVS for metallic lines, as e.g. 
the two C II lines in the red wing of Ha. However, H/3 and H 7 
both clearly show weak variability in the wings. Evaluating 
the variability in Ha is more difficult due to the large num¬ 
ber of telluric lines in this region of the spectrum. In order 
to minimize their influence, only relatively nn-contaminated 
spectra were used in calculating the Ha TVS, with 9/22 
spectra making the cut. There is no variability in the wings 
of this line exceeding that seen in H/3 or H 7 : indeed, there is 
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H He C N 0 Ne Mg Si P S Cl Co TI Cr Fe 




H He C N 0 Ne Mg Si P S Cl Co Ti Cr Fe 



Figure 9. The variability index /van as a function of atomic 
number (top panel) and as compared to the mean EW of the 
line used to measure /var (middle panel). The dashed line is the 
best-fit power law. The bottom panel shows /var.Ci the variability 
index corrected by the hyperbolic fit. 


no evidence of any variation at all in the wings. In order to 
check that discarding spectra was not hiding variability, the 
same observations were used to create test TVS for 'Rj3 and 
Hy: in both cases, the TVS showed variability at a similar 
level as seen in TVS created using the full dataset. 

In order to compare the variability in different elemental 
species. Fig. [9] shows the variability index 


/var — 


(EW) 


( 2 ) 


for all elemental species for which relatively strong, un¬ 
blended lines are available. Measurement of the quantity /var 
was based S^'^-order sinusoids fit to the data (chosen as the 
minimum order necessary to achieve a good fit to the more 
complex EW curves), with uncertainties in the EW and the 
mean EW coming from the uncertainty in the fitting param¬ 
eters. By thus making use of the full dataset, rather than 
simply the two most extreme measurements, uncertainties 
are minimized. 

Weaker lines are expected to be more strongly affected 
by smaller relative changes than will strong lines, and in¬ 
deed, when /var is compared to the mean EW (middle panel 



of Fig.EI, there is an obvious trend of decreasing /var with 
increasing {EW). In order to correct for this trend, we fit a 
purely empirical two-part power law C to the lower bound¬ 
ary of /var, with C = 0.0008(i?W)“^ -I- 0.01(i?W)“°'^. This 
is indicated by the dashed line in Fig. |9] It is plain that the 
differences in variability amongst the lines cannot be due to 
line strength alone. 

The bottom panel of Fig. |9] shows /var,c = /var/C. El¬ 
ements with /var,c = 1 are of course still variable. Eor in¬ 
stance, P disappears almost entirely around phase 0.5 (see 
the dynamic spectra and EWs in Eig. IA2II . and its ex¬ 
tremely high /var reflects this; however, when corrected for 
line strength, it is no more variable than most other lines. 
/var,c thus serves to highlight elements showing particularly 
large line strength variations. In particular. He, C, O, Si, 
and S all possess /var,c > 1- 


6 ABUNDANCES 

We performed an abundance analysis of both HR 2949 
and H R 2948 using the spectrum synthesi s program zee- 
man (iLandstreetl Il988l : IWade et al.l l200lh . zeeman per¬ 
forms polarized radiative transfer under the assumption of 
LTE. Atomic data was taken from the Vienna Atomic Line 
Database VALD2 jPiskunov et akll 19951 : iRvabchikova et ^ 
ll997l:lKuDka et aLlfl999l. 120001) . using an ‘extract stellar’ re- 
quest . Model atmospheres computed with ATLAS9 (iKurucd 
I 1993 I) were used, assuming solar abundances. Surface chem¬ 
ical abundances, as well as usini and TeS, were deter¬ 
mined by directly fitting synthetic spectra to observed 
spectra, using the minimization algorithm described by 
iFolsom et al. l|2012l ). 

The observed spectra for HR 2949 and HR 2948 used 
here were created by averaging all available ESPaDOnS data 
for each star. The averaged spectra were normalized by fit¬ 
ting a low order polynomial through continuum points of 
each spectral order, then dividing the order by that polyno¬ 
mial. This averaging process has the advantage of minimiz¬ 
ing the impact of surface chemical spots, providing results 
closer to the mean global surface chemical abundances. 

Eor HR 2949 the fitting process was performed on 9 
independent spectral windows, covering most of the opti¬ 
cal range (391.5-394, 440-460, 460-480, 490-510, 510-530, 
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Table 4. Absol ute chemical abund ances. Solar abundances 
[X/H ]0 are from lAsplund et alJ l|2Q09lb The number of spectral 
windows used in the analysis is given by N. Abundances rela¬ 
tive to solar are shown in Fig. 1101 For HR 2949, the variability 
quantity /var is also provided (see Fig. 0 . 


El. 

[X/H]q 

HR 2949 

[X/H] /var 

N 

HR 2948 
[X/H] 

N 

He 

-1.07±0.01 

-1.47±0.20 

0.25±0.04 

3 

-1.22±0.30 

1 

c 

-3.57±0.05 

-3.72±0.23 

0.31±0.04 

8 

- 

- 

N 

-4.17±0.05 

-3.41±0.18 

0.27±0.05 

8 

- 

- 

o 

-3.31±0.05 

-3.43±0.33 

0.44±0.08 

8 

-3.12±0.15 

1 

Ne 

-4.07±0.05 

-3.43±0.20 

0.20±0.04 

2 

-3.75±0.40 

2 

Mg 

-4.40±0.04 

-4.11±0.22 

0.07±0.02 

8 

-4.45±0.12 

5 

Al 

-5.55±0.03 

-6.18±0.24 

- 

8 

-5.83±0.30 

1 

Si 

-4.49F0.04 

-4.12±0.28 

0.21±0.05 

8 

-4.55±0.14 

5 

P 

-6.59±0.03 

-6.03±0.30 

2.01±0.83 

1 

- 

- 

S 

-4.88±0.03 

-5.29±0.25 

0.37±0.03 

8 

-4.91±0.30 

3 

Cl 

-6.50±0.30 

-4.94±0.21 

0.34±0.11 

3 

- 

- 

Ca 

-5.66±0.04 

-5.89±0.30 

0.12±0.03 

1 

- 

- 

Ti 

-7.05±0.05 

-5.52±0.30 

0.45±0.14 

2 

-6.81±0.30 

1 

Cr 

-6.36±0.04 

-5.05±0.27 

0.64±0.20 

8 

-6.25±0.30 

2 

Fe 

-4.50±0.04 

-3.83±0.16 

0.21±0.07 

8 

-4.55±0.26 

5 

Ni 

-5.78±0.04 

-5.28±0.19 

- 

8 

- 

- 

Ba 

-9.82±0.09 

-7.79±0.40 

- 

1 

- 

- 

Ce 

-10.42±0.04 

-6.64±0.20 

- 

1 

- 

- 

Pr 

-11.28±0.04 

-8.16±0.21 

- 

8 

- 

- 

Nd 

-10.58±0.04 

-7.43±0.25 

- 

8 

- 

- 


530-550, 550-570, 570-600, and 600-650 nm), excluding H 
lines, the strongest He lines, and regions contaminated by 
telluric features. The results were then averaged to produce 
the final best-fit values, with the standard deviation taken 
as the uncertainty. For elements with constraints from less 
than 4 spectral windows, the uncertainties were increased 
to include the full window to window scatter, potential nor¬ 
malization errors, and the impact of any potential blending 
lines. We used a hxed logp = 4.1, based on the Balmer line 
analysis. 

While a major strength of zeeman is that it can include 
the effect of magnetic fields and perform polarized radiative 
transfer (integrated across a stellar disk), doin g so increases 
the r un time by 3 orders of magnitude (see I Folsom et al.l 
l2012h . Due to computational limitations we can either per¬ 
form an extremely accurate analysis of a few elements and 
spectral lines, or a more approximate analysis of many more 
elements and spectral lines. Since we are interested in diag¬ 
nosing the full range of peculiarities in HR 2949, as well as 
determining Teff based on lines with a wide range of ioniza¬ 
tion states and excitation potentials, breadth is more impor¬ 
tant for this analysis than extreme precision. Therefore the 
analysis of HR 2949 was performed using a non-magnetic 
model, with broadening due to Zeeman splitting approxi¬ 
mated by 2 kms“^ of microturbulence. In order to check 
the validity of this approximation, we conducted a second 
analysis in the 490-510 nm spectral window for He, N, Si, S, 
Cl, Cr, Fe, Ni, and Nd, assuming Ba = 2.7 kG in the geom¬ 
etry derived in Sect 7.2 for phase 0 (and using the same fit¬ 
ting procedure as above). The abundances determined using 
magnetic vs. non-magnetic models typically agree to better 
than 0.02 dex, and always better than 0.05 dex. Thus, the 
impact of the magnetic field is much smaller than our other 
uncertainties. 

We find HR 2949 has a best fit Tefi of 18.3 ± 0.8 kK, 
and a usini of 63 ± 3 kms“^. This Teg is in agreement 
with our photometric and ion strength ratio measurements 



Figure 11. The corrected variability index /var.C as a function 
of relative solar abundance. Vertical and horizontal dashed lines 
indicate, respectively, solar abundance, and the mean level of vari¬ 
ability for a line of the strength used to measure /var. 

in Section 4. We find strong photospheric chemical peculiar¬ 
ities in HR 2949, reported in Table HI and plot ted relative 
to the solar abundances of lAsplund et al] (l2009h in Fig. 1101 
We find strong overabundances of Fe-peak elements, by 1- 
2 dex above solar, and very strong overabundances of the 
rare earths Ce, Pr and Nd, by 3-4 dex above solar. Con¬ 
versely, we find an underabundance of He, as well as Al and 
S. We also hnd strong overabundances of Ba, Cl, and pos¬ 
sible overabundances of Ne and P. This pattern of chemical 
abundances is characteristic of magn etic Ap/Bp stars (e.g. 
iFolsom et al.l[2007l : lBailev et al.ll20ld l. thus we conclude HR 
2949 is chemically a Bp star. 

We performed a similar abundance analysis on HR 
2948, for comparison. Somewhat larger spectral windows 
were used (440-480, 490-530, 530-570, 570-610, and 610- 
650 nm), due to the higher usini producing weaker, more 
blended lines. A logg = 4.24 was used, based on the Balmer 
line analysis. We find a TeS of 13.6 ± 1.2 kK and a usini of 
117 ± 2 kms”*^. We attempted to fit microturbulence as 
well, however the value is poorly constrained, since the large 
V sin i leads to a lack of unblended weak lines (on the linear 
part of the curve of growth). Thus we only constrain micro¬ 
turbulence to be ^ 1.2 kms“^. The chemical abundances 
we find for HR 2948 are entirely consistent with solar (see 
Table H] and Fig. [TO} . Therefore we conclude that HR 2948 
is a chemically normal B star. 

A full comparison of variability and abundance cannot 
be performed, as strong, unblended lines do not exist for 
many of the elements examined in the abundance analy¬ 
sis, particularly the rare-earth elements. For those elements 
which can be so analyzed. Fig. Illl compares the abundances 
to the quantity /var.c from equation [2] and Fig. [9] There is no 
clear relationship between the variability of a given element 
and its abundance. The spectral lines of elements with so¬ 
lar, sub-solar and super-solar abundances are in some cases 
highly variable, and in other cases not. 
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Figure 12. Representative LSD profiles for HR 2948 (top) and 
HR 2949 (bottom). The HR 2949 LSD profile was made using 
a line mask including all metallic lines and the 2010-01-01 ES- 
PaDOnS observation. Stokes I (bottom) is in black, Stokes V 
(top) in red, and the diagnostic null N (middle) in blue. Inte¬ 
gration ranges used to measure {Bz) are indicated with dotted 
vertical lines. Note the strong Stokes V signature in HR 2949, in 
contrast to that of HR 2948, which is consistent with noise. 


7 MAGNETIC FIELD 


In order to maximize the signal-to-noise ratio, the ES- 
PaDOnS observations were analyzed usin g Least Squares 
Deconvolution (LSD) (|Donati et al.lll997|) as implemente d 
in the iLSD package developed by Kochukhov et al. 

As in the ab undance analysi s, line l is ts were downloaded 
from yALD2 JPisk unov et all Il995l : iRvabchikova et Tl 
I 1997 I : iKupka et al.l I 2 OOC 1 I I. For HR 2949, the line list was 
constructed for a star with Tefj^ 18 kK, log(;= 4.0, and a 
line-strength threshold of 10% beneath the continuum, yield¬ 
ing 462 lines. The line list for HR 2948 was constructed for a 
star with Tefr= 14 kK and logg= 4.0, and a 10% threshold, 
yielding 561 lines. 

The line lists were cleaned by removing all H and He 
lines as well as all lines blended in their wings, along with 
all lines strongly contaminated with telluric lines. While in 
most cases the effect on Stokes V of these blends is negligi¬ 
ble, cleaning the lists in this fashion ensures that the mean 
Stokes I profiles are not affected by lines with strong depar¬ 
tures from the mean photospheric line profile. The remaining 
lines were then empirically adjusted such that their depths 
reflect the observed line intensities. 


Ha 



Phose (P = 1.90871 d) 


Figure 13. Single-line longitudinal magnetic field measurements 
using Ha (solid black circles). The solid (black) curve is the best- 
fit sinusoid; the dashed (red) curve is the best-fit offset dipole 
model with a = 0.33. 


For HR 2949, LSD profiles were created using a mask in¬ 
cluding all metallic lines, as well as single-element line masks 
for C, N, O, Ne, S, Mg, Si, and Fe. A He line mask was also 
tested. No single-element masks were created for HR 2948, 
although He was excluded in order to avoid distortion of the 
Stokes I profile. The number of lines in each mask is given 
with the longitudinal field measurements in Table IBTI Rep¬ 
resentative LSD profiles created using metallic masks are 
shown for both stars in Fig. 1121 

To determine whether or not a magnetic field is de¬ 
tected, a statistical test is applied JPonati et al.l[l99^ . A 
detection is considered ‘definite’ if the signal within the 
Stokes V line profile departs from the noise, as diagnosed by 
both Stokes V outside the line profile, and by the diagnostic 
null N, with the False Alarm Probability FAP < 0.001%. 
If FAP > 0.1%, the observation is formally considered a 
non-detection. Values between these limits are considered 
marginal. Integration ranges, determined by the extent of 
the absorption lines, are indicated by dotted vertical lines 
in Fig. 1121 For HR 2949, all measurements are definite de¬ 
tections, with FAP <C 0.001. For HR 2948, all individual 
observations are non-detections, as is the co-added observa¬ 
tion. 

The longitudinal magnetic field (Rz) was measured 
from the LSD profiles b y taking the fi rst-order moment of 
the Stokes V profile (e.g. iMathvjl 19891 ): 


(B^) = -2.14 X 10“ 


f vV{v)dv 

Xogoc f [7c - I{v)] dv’ 


(3) 


where v is the Doppler velocity in kms“^, and Ao and go 
are the normalization values of the wavelength and Lande 
factor used to scale the Stokes V profile. For all LSD profiles 
the same parameters were used to extract the LSD profiles 
and to measure (Rz): Ao = 500 nm and go = 1.2. The in¬ 
tegration ranges are shown as dotted lines in Fig. 1121 The 
same measurement can be applied to the diagnostic N pro¬ 
file, yielding an equivalent ‘null longitudinal magnetic field’ 
(Nz), with which (Rz) can be compared. 

LSD profile measurements of (Rz) for HR 2949 are given 
in the Appendix in Table IBTI For HR 2948, the single mea- 
















































Magnetic field and speetral variability of HR 2949 11 


surement gives {Bz)= —19 ± 35 G and {Nz)= —69 ± 35 G, 
consistent with a non-detection. 

Because Zeeman signatures are visible in many of HR 
2949’s spectral lines, single-line measurements of (Bz) were 
also performed for a selection of lines: Ga ll 393.4 nm and Ti 
II 457.2 nm (since there are too few lines of these species to 
perform LSD); Si ill 455.2 nm and Fe ll 516.9 nm (for com¬ 
parison of single-line results to LSD results), and Hq (since 
the relatively uniform distribution of this element across the 
photosphere should result in Stokes V profiles that are most 
representative of the real global magnetic field geometry). 
These measurements, together with the rest wavelengths Ao 
and Lande factors g, are given in the Appendix in Table [B^ 

The Hq (Bz) measurements are shown phased with the 
rotation period in Fig. 1131 The best-fit sinusoid is overplot¬ 
ted as a thick black line. The overall shape of the (Bz) curve 
is somewhat non-sinusoidal, being sharply peaked near the 
positive pole and relatively broad near the negative pole. 
However, only two points are incompatible with a sinusoidal 
fit at the Sa level. 


7.1 Magnetometry using different elements 

The LSD Stokes I and V profiles of HR 2949 created from 
He, metal, O, Si, S, and Fe masks are shown in Fig. 1141 The 
LSD profiles are stacked vertically, with the vertical offset 
proportional to the rotational phase. The LSD profiles cre¬ 
ated using all metallic elements are most similar to those 
created with Si and Fe masks, suggesting the former is dom¬ 
inated by these lines. 

Comparing the Si and Fe LSD profiles, Fe lines are con¬ 
sistently more complex than those of Si lines. This is likely 
a reflection of the greater complexity of the photospheric Fe 
abundance distribution, as is apparent in both the Stokes I 
profiles and the dynamic spectra in Figs. FtI and IA2I 

The strongest difference between the Fe and Si lines 
(right column of Fig. 1141) . and the He, O, and S lines (left 
column of Fig. [Mil, occurs near the positive magnetic pole. 
Here, the Stokes V signatures of the latter group disappear 
almost entirely, while those of the former group show cross¬ 
over signatures followed by a reversed polarity. 

O and S are, however, quite different from He: while He 
shows a simple S-shaped Stokes V profile near the negative 
magnetic pole, O and S both display more complex Stokes 
V profiles. 

Fig. [T^ shows mosaics of phased LSD (Bz) measure¬ 
ments from Table IBII Best-fit sinusoids are shown in red. 
Measurements of He, O, and Ne are apparently well-fit by a 
sinusoid; however, they do not reverse polarity, while most 
other elements (except Mg) do. Of the remaining elements, 
some, such as C, Si, and Fe, are clearly non-sinusoidal. 

Fig. lfbl shows a mosaic of (Bz) measurements performed 
on individual spectral lines (see also Table [B^ . The Si ill 
and Fe ll measurements are similar to the LSD measure¬ 
ments for these species, albeit with greater scatter and larger 
error bars. While the best-fit sinusoid to the Ti ll 457.2 
nm measurements indicates polarity reversal, the error bars 
of most measurements near the positive pole are consistent 
with zero. Ga ll, however, clearly belongs to the non-polarity 
reversing group. 

Comparing spectral and magnetic diagnostics, there ap¬ 
pears to be an absorption excess of He associated with 
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Figure 14. Phased LSD profiles, with the continuum level set to 
the rotational phase. The Stokes V scaling factor is relative to 
Stokes I. 


the longitude of the negative magnetic pole, and a weak¬ 
ening of absorption associated with the positive pole. This 
leads to an apparent weakeni ng of the positive ma gnetic 
pole when viewed in He lines. lYakunin et al.l ll2014h have 
reported similar behaviour for HD 184927, and have shown 
that by self-consistently mapping the abundance distribu- 
tions and magnetic field using Zeeman Doppler Imaging (e.g. 
IPiskunov fc Kochukhovll2002l ~) the divergent (Bz) measure¬ 
ments from various chemical species can be reproduced with 
a single magnetic model. 
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Figure 15. LSD profile longitudinal magnetic field measurements 
for single-element masks. 
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Figure 16. Single line longitudinal magnetic field measurements. 


7.2 Magnetic Field Geometry 

IPrestonl lll967ll showed that, for a dipolar magnetic field in¬ 
clined at an angle P from the rotational axis and an inclina¬ 
tion i of the rotational axis from the line of sight, i and /3 
are related by 


tan P 



where r is the ratio 

|Bo| 

\Bo\+BS 


(4) 

(5) 


and Bq and Bi are dehned by the sinusoidal fit to the phased 
(f?z) variation: 


Bz = Bo -\- Bi sin 27r((^ — po)- (6) 

As the (Bz) variation differs depending on the spec¬ 
tral lines with which it is measured, so too do Bo and Bi. 
We base our primary analysis upon Hq measurements, as 
these are expected to be the least affected by LPV. Fit¬ 


Table 5. Sinusoidal fitting parameters Bq and Bi, and resultant 
dipolar magnetic field model parameters Bo and 0 for single¬ 
element (Bz) measurements (see eons. [BHTll. 


Element 

Bo 

Bi 

Bd 

P 


(G) 

(G) 

(kG) 

(°) 

H 

-278±14 

376±18 

2-4to.2 

50±16 

He 

-514±14 

291±19 

q o+J.-'i 
'^•^-0.4 

27±16 

C 

-396±47 

716±60 

^•^-0.3 

57±15 

N 

-144±64 

1037±82 

^•-*^-1.0 

80±5 

O 

-450±20 

-69±26 

^•^-0.7 

8±6 

Ne 

-533±27 

193±36 

q 

'^•^-0.6 

19±13 

Mg 

-450±37 

-454±50 

o 4+u.» 

42±17 

Si 

-271±25 

726±32 

o q+1-2 

66±11 

S 

-289±42 

575±55 

q q-|-U.C5 

'^•'^-0.3 

59±14 

Ca 

-495±57 

506±73 

q Q+u.y 

42±17 

Fe 

-209±29 

632±38 

o o + l-l 
'^•‘^-0.4 

68±10 


ting parameters for different single-element measurements 
are tabulated in Tabled 

The best-fit sinusoid to the Ha measurements has 
Bo = - 278 ± 14 G and Bi = 376 ± 18 G, giving 

r = — 0.15 ± 0.05. With the inclination determined 

above, we then have /3 = 50 ± 16°. 

The strength of the magneti c dipole can then be derived 
from the relation (|Prestonlll967l i 

Bd = B™ ^ 20(^3~*-^) Pcosi + sin P sin , (7) 

where u is the limb darkening coefficient, which we take to 
be 0.32, the mean of the values obta ined from the tables 
calcul ated by Ivan Hammd (Il993li and iDiaz-Gordoves et al.l 
l|l995li . This yields Ba = 2.4+°;3 kG. 

If metal lines are used instead, equations [ 6 ] and □ of 
course give very different results for P and Ba, which dif¬ 
fer outside the formal uncertainties. These are summarized 
in Table [5l This demonstrates that, in the presence of sur¬ 
face abundance inhomogeneities, the use of particular sets 
of spectral lines can introduce significant systematic error 
into the determination of the magnetic field strengths and 
geometries of Bp stars. 

Fig. ll7l shows a comparison between synthetic Stokes V 
and observed Ha Stokes V profiles. The model synthetic pro¬ 
files were calculated using a disk integrati on method similar 
to that described bv IPetit fc Wada (l2012li . We examine two 
models: a dipolar model (solid blue lines) using the param¬ 
eters above, and the best-fit offset dipole model (dashed red 
lines). The observed profiles (black points) exhibit a much 
simpler variation than is seen in He or metallic lines (com¬ 
pare to Fig. El). The dipolar profiles are able to reproduce 
the observed profiles quite well around the negative mag¬ 
netic pole (rotational phases 0.35-0.5), however the fit is not 
as successful in the vicinity of the positive pole (rotational 
phases 0 . 8 - 1 . 2 ). 

To determine the parameters of the offset dipole model, 
we employed a goodness-of-fit test on the Ha (Bz) mea¬ 
surements. The (Bz) curve is shown in Fig. El Th® model 
parameters are P = 60°, Bd = 2.4 kG, with the dipole offset 
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Figure 17. Stokes V profiles from Ha (black) arranged in order 
of phase. Solid blue lines are synthetic Stokes V spectra created 
for HR 2949’s stellar parameters and a dipolar model matching 
the best fit model to the Ha (Bz) measurements. The model 
agrees well with observations near the negative magnetic pole 
(phase 0.5), but cannnot reproduce observations near the positive 
magnetic pole (phase 0.0). Dashed red lines are synthetic profiles 
created for an offset dipole model with a = 0.33: while this model 
provides a somewhat better fit near the positive pole, the fit near 
the negative pole is substantially worse than the dipolar model. 


by a = 0.33R* along the magnetic axis. As is clear from Fig. 
El even this large offset is unable to satisfactorily reproduce 
Stokes V: while the fit is somewhat improved near the pos¬ 
itive pole, it is somewhat worse near the negative pole. In 
both cases, the fit is especially poor between phases 0.6 and 
0.9. 

While H lines are not as subject to differential flux re¬ 
distribution as He and metallic lines, it is clear from both dy¬ 
namic spectra and EW measurements that H is so affected, 
likely due to variations in the local H to He ratio (see Fig. 
0 . Furthermore, it is precisely at those phases that are most 
poorly fit that the line strength is at its strongest. From this 
we conclude that H line Stokes V profiles are likely subject 
to warping by differential line strengths, and that Magnetic 
Doppler Imaging (MDI) will be necessary to reconstruct the 
magnetic field topology. 

8 DISCUSSION 
8.1 Binarity 

There is a discrepancy of -^100 Myr in the apparent ages of 
HR 2948 and HR 2949 (see Table [3]). Since binary pairs are 
expected to be coeval, this motivates a re-examination of 
the assumption that the stars are a binary pair, as, if this is 
not correct, some of the stellar parameters found above may 
need to be revised. Here we explore four different scenarios 
to attempt to explain this age discrepancy: 

1) The Hipparcos parallax is unreliable, and hence the dis¬ 
tance and luminosities found above are incorrect. 

2) The isochrones used to determine the ages are inappro¬ 
priate: pre-main sequence (PMS) isochrones should be used 
instead. 

3) The stars are not gravitationally bound: their proximity 
on the sky is a chance coincidence, and they are not in fact 
at identical distances. 

4) The stars are not coeval, but form a binary pair due to 
gravitational capture. 

Keeping Teg the same and setting the minimum lumi¬ 
nosity of HR 2949 at the ZAMS, the minimum distance to 
the star is 112±6 pc. If HR 2948’s luminosity is recalculated 
at this distance, it’s age is 89jl3j Myr, i.e. the star is well off 
the ZAMS, and the age discrepancy is not resolved. Clearly, 
if instead the distance to the stars is greater than that deter¬ 
mined from the Hipparcos parallax, the discrepancy is im¬ 
possible to resolve: HR 2949 would have to be much cooler 
to be on the same isochrone as HR 2948. Therefore scenario 
1 ) does not resolve the problem. 

To e xplore scenario 2), we use the PMS isochrones pre¬ 
sented bv lTognelli et al] (l201l| j. Only the 1 Myr isochrone is 
consistent with the luminosities and effective temperatures 
of both stars. This isochrone is shown by a dash-dotted line 
in Fig. |4l In this scenario, HR 2949 is essentially at the 
ZAMS, while HR 2948 has not yet reached it. 

However, HR 2948/9 is not located in or near to any 
star-forming regions, as might be expected if the system is 
indeed 1 Myr old. Furthermore, HR 2948 has no Ha emis¬ 
sion. Therefore it seems unlikely that the system is newly- 
formed. 

To explore scenario 3), we note first that the radial ve¬ 
locities are identical within error bars, as expected for a 
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weakly gravitationally bound system. The proper motions of 
the two stars differ by 8.8±2.6 mas yr ~^ in right ascen sion 
and 4.8±2.6 mas yr~^ in declination (iHog et al.ll200fll : see 
also Table [3]), for a total relative proper motion of 10.0±3.5 
mas yr“^. At the Hipparcos distance of 139 pc, this cor¬ 
responds to a relative motion of b.blj'g kms“^. The stars 
are separated on the sky by 7.3”, corresponding to a phys¬ 
ical separation of lOldljg® AU. With the masses above, 
and assu ming a circular orbit, the orbital velocity is then 
Worb ~ G(Mi -I- M 2 )/r = 3.0±0.3 kms“^. In order for the 
orbital velocity to equal the relative motions, the semi-major 
axis of the system would have to be a = 90itgg° AU, clearly 
inconsistent with the observed separation. The proper mo¬ 
tions thus argue against the two stars being gravitation¬ 
ally bound. This analysis assumes the orbital plane is in the 
plane of the sky. If the inclination of the orbital axis from 
the line of sight is substantially greater than 0°, the physical 
separation would of necessity be larger, implying a lower Ke- 
plerian velocity. At the same time, the deprojected relative 
velocity would be larger, thus increasing the discrepancy. 
However, the disagreement between the orbital velocity and 
the proper motions is relatively small, so we do not view this 
as conclusive evidence that the two stars are not associated. 

Regarding scenario 4), it is worth noting that the system 
is a wide binary, and not strongly gravitationally bound, 
making a gravitational capture scenario more plausible. 


8.2 Variability 

The most likely explanation for the LPV seen in HR 2949’s 
He and metallic lines are surface abundance variations. 
Chemical spots are commonly see n in Ap/Bp stars (e.g. 
lMichau3ll97nl : iMichaud et al.lll98ll 'l , indeed the presence of 
such spots is one of the most reliable indirect indicators of 
magnetism amongst early-type stars. Differences in the cor¬ 
rected variability index /var,c are almost certainly related 
to this: elements with /var.c > 1 likely show more extreme 
surface abundance variations than elements with /var,c = 1- 

Comparing Figs. [T] [TS] and 1151 the strongest He, O, 
and S line profiles are associated with the negative magnetic 
pole, while the positive magnetic pole is associated with ab¬ 
sorption minima. The weakness of these lines near the pos¬ 
itive pole may explain why their Zeeman signatures seem 
to disappear almost entirely near the positive pole (Figs. 
1141 and IT^ . The intensity of the regions of the photosphere 
associated with the positive pole is reduced, decreasing the 
contribution of the magnetic pole to the polarized flux. 

A similar effect in metallic lines is likely the reason for 
the observed inconsistency in (Hz) as measured using differ¬ 
ent elements. The strong asymmetry observed in the LSD 
Stokes V profiles (see Fig. I14II is associated with asymme¬ 
tries in Stokes I. It is also possible that the apparent de¬ 
partures of (Hz) from the simple sinusoid expected for a 
purely dipolar magnetic field are likewise a consequence of 
this effect, with greater relative weight given to the positive 
magnetic pole magnifying ( Hz) at these phases. S uch a n ef¬ 
fect was first suggested by iBorra fc Landstreetl lll977li for 
the magnetic Ap star CVn, who noted that magnetom- 
etry conducted using H lines did not show the non-dipolar 
character of the (Hz) measurements conducted using metal¬ 
lic lines. 

Conversely, the departures of Ha measurements from 


purely sinusoidal behaviour, along with the failure of a dipo¬ 
lar model to reproduce the Stokes V profiles near the posi¬ 
tive magnetic pole in this line, point towards a non-dipolar 
component to the magnetic field. 

Examination of the Hp magnitudes in Fig. [T] show 
that the star is brightest between phases 0.0 and 0.2, and 
dimmest between phases 0.4 and 0.5. This corresponds quite 
well to the line-strength extrema of He and Si. It thus seems 
likely that the photometric variability of HR 2949 is a con¬ 
sequence of the variable brightness of photospheric abun- 
danc e patches, typical of magnetic chemically peculiar stars 
(e.g. IWraight et al.l[2012|j . Efforts to model the photomet¬ 
ric variations of some Bp stars by taking into account the 
variable optical brightness caused by flux redistribution due 
to chemical spots have successfully accounted for a majority 
or all of the light curve variabilit y in many cases, e. g. for 
the He-strong B p star HD 37776 (iKrticka et aDl2007l L the 
Ap star t UMa Ishulvak et al.l[2010bl L the Si star HR 7224 
jKrticka et al.ll2009l L and the Ap star CU Vir llKrticka et al.l 

I 2 OI 2 II . 


While the strongest H variability is seen in the core 
of the lines, and is almost certainly also a consequence of 
photospheric abundance variations, the origin of the LPV 
observed in the wings of the H Balmer lines (see Figs.Qand 
EJ is not so obvious. 

Early-type Bp stars often show emission in the 
wings of their Balmer lines origina ting in their corotat¬ 
ing ‘centrifugal m agnetospheres’ (e.g. iTownsend et al.|l2005l : 
IPetit et al.l I 2 OI 3 I I . The magnetospheric parameters of HR 
2949 are derived in the following subsection. Here, we 
note that the Balmer line variability is extremely weak; 
shows no signs of the blue-to-re d variation associated 
with magnetospher ic plasma (e.g., iGrunhut et al.l l2012cl . 
iRivinius et al. 2013lL and shows no sign of a Balmer decre¬ 
ment I Williams fc Shipman|[l988h . as would be expected for 
a transluscent circumstellar plasma (indeed, Ha is less vari¬ 
able than the higher-numbered Balmer lines, precisely the 
opposite of expectations for circumstellar emission). There¬ 
fore we conclude that the Balmer wing LPV is not a conse¬ 
quence of a magnetosphere. 

Another possibility is that variation in the wings is re¬ 
lated to the chemical spots. This could arise due to different 
degrees of line blanketi ng at different points of the surface 
dKhan fc Shulvakll2007h . producing small differences in at¬ 
mospheric structure in different regions of the photosphere. 
In this case, the Balmer wing s trength should correlate with 
Fe, and anticorrelate with Si (IKhan fc ShulvaH l2007l L ap¬ 
proximately matching observations. 

A third possibility related to chemical spots is that the 
LPV is a pressure effect originating in local abundance en¬ 
hancements or depletions. Such an eff ect is suspected in th e 
magnetic Bp star HD 184927 (B2V. [Yakunin et aLllioi^ L 
The variations in HR 2949’s wings are much weaker than 
those in HD 184927, however HD 184927 is a He-strong star, 
with large surface abundance variations. 

The final possibility is that the LPV is a consequence of 
variable magnetic pressure originating due to induced elec¬ 
tric currents in the photosphere. This phenomenon has been 
suggested as the reason for Balmer li ne LPV in the He-wea k 
Si star HR 7224 (B9p, Hd< 0.4 kG iLehmann et al.ll2007l L 
In this case chemical spots were specifically excluded as the 
possible source of the LPV. Model atmospheres including 
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Table 6. Mag netospheric param eters calculat e d ass uming mass- 
loss rates from IVink et alJ <1200 ih and iKrtickal ll2014l 'l . 


Rk (R*) 

3.5 

+1.7 

-1.2 


W 

0.14 

-HO.10 
-0.07 


Mass-loss prescription 

Vink 


Krticka 

logM (M© yr-l) 



-lO-Stai 

Uoo (kms“l) 

1070t“ 



log 

4.7t°;« 


5-5la? 

Ra (R*) 




log (Ra/Rk) 

0.63±0.35 


0-82l°« 

TM (Myr) 

S-Oti g X lO'^ 

1 .8; 

tag X 10® 

Tj (Myr) 

11±1 


32±13 

(Myr) 

22±3 


63±35 


the Lorentz force were shown to provide the b est explanation 
for B aimer line LPV in the AOp star 6 Aur llShulya k et alJ 
l200it l and the B6p star 56 Ari JShulvak et al.l2010al i. both of 
which host moderately strong photospheric magnetic fields 
{Bd^ 1 kG). HR 2949’s Balmer line wing LPV is similar 
in character to that seen in HR 7224, 9 Aur, and 56 Ari, 
suggesting that a similar mechanism may be at work. 


8.3 Magnetosphere 


When contained within a magnetic field, the spherical 
symmetry of a stellar wind is strongly modihed if the 
energy density of the magneti c field exceeds the kinetic 
energy density of the w ind llBabel fc Montmerld 1 19971 : 
lud-Doula fc Owockil l2002ll . Confined plasma is forced to 
flow along closed magnetic loops. At the tops of the 
loops, the plasm a collides, releasing high-energy X-rays in 
wind s hocks (e.g. Bal^fc Montmerlelll99'4 : lud-Doula et al.l 


I 2 OI 4 I : iNaze et'S! l201^ The plasma may also become 
detectable as rotationally-modulated emission in wind- 
sensi tive spectral lines, particularly in the ultraviolet 
fe.g. IShore et al. 1990j: Smith fc Grootell200l] : ISchnerr et al.l 
I 2 OO 8 I: lHenri~ .eL.alJ ( 201 ^) , but also in optical lines such 
as Hg (e.g. WalbornI |l 974 iPed ersenI Il 979l: iLeone et al 

iill 


2013). 


20K: [&unhut et al.ll2012al : Rivinius et ahl 2 OI 3 O Petit et al 


As discussed above, we do not detect any sign of 
Ha emission in HR 2949’s spectra. 

There have been no published X-ray observations of HR 
2949. 


There is a single high-dispersion lUE spectrum, ac¬ 
quired with the large aperture, a 10x20” oval which almost 
certainly included HR 2948. The flux is much too high to 
be associated with a single star of V ~ 4.5. However, it is 
compatible with the combined magnitude of the binary sys¬ 
tem. Using the ephemeris in equation[T] the lUE observation 
was acquired at a rotational phase of 0.27±0.12. Thus, the 
observation was either acquired during cross-over, at which 
point the magnetic equator is visible, or close to the neg¬ 
ative magnetic pole. Approximating the magnetosphere as 
a torus in the plane of the magnetic equator, in the former 
case it should be detectable as enhanced absorption in the 
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Figure 18. The single lUE high-dispersion observation (black), 
with photospheric model spectra for HR 2948 (dashed blue), HR 
2949 (dotted magenta), and the combined system (solid red). The 
three panels show wind-sensitive UV doublets; short vertical lines 
indicate the doublets’ rest wavelengths. 


cores of wind-sensitive resonance lines, as the plasma will 
be occulting the disk. In the latter case it should appear as 
emission peaking at high velocity, as the plasma i s projected 
on either side of the disk llSmith fc Groot jl200ll ). 

Fig. [18] shows the Si iv, G iv, and Al iii doublets in 
the lUE observation, as compared to a photospheric model 
created for HR 2948-1-2949 by combining zeeman/atlas9 
models generated using the stellar parameters in Table 
[3| line lists from VALD, and the abundances in Table |4| 
Th is version o f zeem an included modifications introduced 
by iLandstreed (l201ll ). who added extra continuous opacity 
sources for the UV. The models were rotationally broadened 
using the usini of each star, and convolved to the FWHM 
of the lUE data (0.1 A). Fluxes were weighted by stellar 
radii (4.2 77© and 3.7 7?© for HR 2948 and HR 2949, respec¬ 
tively), and dereddened assuming the same small reddening 
{E{B — V) = 0.04) used for the photometric analysis in 


Section 4. 

These three doublets are sensitive to the stellar wind, 
and are commo nly seen in emiss i on in magnetic early- 
type stars (e.g. . Barker et al.l 1983: Smith fc Groo^l200ll : 
lOskinova et al.l l201ll : iHenrichs et ahl l2013l ). In this case, 
there is no evidence of emission in any of these lines. 

There is thus no direct evidence for the existence of 
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a magnetosphere. However, its expected properties can be 
investigated based upon HR 2949’s magnetic and stellar pa¬ 
rameters. The magnetospheric parameters derived below are 
summarized in Table [6l 


8.3.1 Theoretical Characteristics 

The theo ry of magnetic wind confi n ement is treated in 
detail by iBabel fc Montmerld (Il997l h Ind-Doula fc Owockil 
(|2nn2h . and lud-Doula et al.l ll2008h . The specific case of 
the rigidly rotating or ‘centrifugal magnetosphere’ found in 
magnetic mid B-typ e stars snch as HR 2949, is given by 
iTownsend fc Owockil (l2005ll . Summ aries of magnetic wind 
confinement theory are provided bv IPetit et al. ll2013lf and 
lYakunin et al.l (120141 ') , whose basic outline we follow in deter¬ 
mining the following magnetospheric parameters: the mag¬ 
netic wind confinement parameter (which gives the ratio 
of magnetic to wind kinetic energy density); the Alfven ra- 
dins Ra (which gives the maximum radius of closed mag¬ 
netic loops); the rotation parameter W (the ratio of the 
equatorial to the critical rotational velocities); the Kepler 
radius Rk (the radius at which centrifugal force balances 
gravity); and the spindown timescale rj. 

To determine rjt we must hrst estimate the wind mo¬ 
mentum Mvoo, which we evaluate from the stellar param- 
eters determined above (Table 0, using both the recipe of 


1.31, with logjQ Ra/Rk likely belo w 1 using 


loss rates and lik ely above 1 using iKrticka mas s-loss rates. 


Vink et ahlmass 


Vink et ^ lll999l . l2000l . booil l , and the mass-loss rates of 
Krtickal (l2014l l. who tailored the calculations specihcally for 


chemical ly peculiar s tars. 

The I Vink et al.l recipe yields log (M) = —9.7jlQ'3 and 
Voo ~ 1070la n kms“^. 

iKrtickal ’s mass-loss rates are highly sensitive to abun¬ 
dances. The elements to which the wind is most sensitive 
(C, N, O, Si, S, and Fe) are either consistent with solar 
abundances (C and O), slightly underabundant (S), or only 
slightly overabundant (N, Si, and Fe). Moreover, while in 
general increased abundances increase Mv^a, at low Teff in¬ 
creased Fe abundances decrease the wind momentum. To 
make the most use of IKrtickal ’s tables, we therefore in¬ 
terpolated their calculated values of M and Uoo accord¬ 
ing to the abundances in Section 6, then took the mean 
of the values calculated via the various tables. These are 
log (M) = —lO.SjlQ ® Mq yr“^ and Voo= 1070jlggg kms“^. 
The uncertainty in M is largely a function of the uncertain¬ 
ties in Teff and R«, as there is little scatter between results 
for different elemental tables. The uncertainty in Vao reflects 
the enhan cement of th is parameter by the Si overabundance. 

Using l^nk et al.l mass-loss yields logr;, = 4.71 q ® » 1, 
while using iKrtickal mass-loss gives log» 7 « = S.Slj'g. In ei¬ 
ther case, ? 7 * ^ 1, therefore the wind is magnetically con¬ 
fined. The Alfven radius is then Ra = ISjigT?* or Ra = 
23jig'll?*, respectively, i.e. in either case Ra > 107?,, and 
HR 2949 is predicted to possess a magnetosphere of signifi¬ 
cant spatial extent. 

Taking the equatorial rotational velocity, radius, and 
mass found in Section 4 yields Rk= Z.S/WR, and W = 
0.141q;(,°, respectively. Since Rk<Ra, HR 2949 is expected 
to possess a centrifugal magnetosphere. 

The ratio logjg Ra/Rk provides a dimensionless 
measure of the magnetosphere’s volume. In this case, 
logj^g 7Za/7?k spans a wide range of values, from 0.28 to 


IPetit et al.l (l2013l l. who base their analysis upon ivink et al.l 
mass-loss rates, give a value near unity. Our lower value is 
due to the downward revision of 73d, and the upward re¬ 
visions of 7? ,, log I/, and Teff , as compared to the values 
reported by IPetit et al.l (l2013l l: the weaker magnetic field, 
and stronger wind, lead to a reduced magnetic confinement 
volume. 

The magnetic wind c onfinement theory developed by 
lud-Doula fc Owockil (|2002| ~I deals with the simplest magnetic 
topology, i.e. a dipole aligned with the rotational axis. One 
might therefore wonder if the substantial obliquity of HR 
2949’s magnetic field with respect to the rotational axis may 
be the reason for the lack of an observable magnetosphere. 
However, the semi-analytic Rigi dly Rotating Magnetosphere 
(RRM) formalism developed bv lTownsend fc Owockil (l2005l i 
is able to predict the circumstellar plasma distributions 
and associated observational diagnostics for arbitrary val¬ 
ues of /3, and has been successful in qualitatively explain¬ 
ing the rota tionally modulated Ha emission seen in 5 Ori 
C (B3 Vp, Leone et al.| I 2 OI 0 I I and HD 176582 (B5 IV, 
iBohlender fc MoninI 201l] l. both of which have /3 close to 
90'’. The obliquity of HR 2949’s magnetic field is therefore 
an unlikely explanation for the lack of an observed magne¬ 
tosphere. 

A major consequence of the coupling between stel¬ 
lar winds and magnetic fields is rapid angular momentum 
loss, as the corotating pla s ma acts as an extended moment 

arm llWeber fc Davislil967l : lud-Doula et al.ll2009l: IPetit et al.l 


l2013li . From the internal structure models of Claret (200J), 
the moment of inertia / 0.06 for a 6 Mq star at an age of 

-^11 Myr, with very little variation during its previous evo¬ 
lution. Taking /, the mass-loss timescale tm, and Ra to be 
constant, this then yields a spindown timescale rj = 11 ± 1 
M yr or 32 ± 13 Myr, depending on whether the lVink et al.l 
or IKrtickal mass-loss rates are used. 

We can use rj and W to infer the st ar’s spindown age ts 
(lud-Doula et al.l [20091 : IPetit et al.ll2013l ~). However, we note 
that the assumpti on that the initial rotational parameter 
Wo = 1.0, made bv IPetit et al.l (I 2 OI 3 II in their estimation of 
the maximum spindown age, is implausible. While magnetic 
torques are unable to spin down massive magnetic YSOs 
(iRosen et al.l l2012l ~). gravitational torques during star for - 
mation enforce an upper limit of Wo = 0.5 (iLin et al.ll201ll l. 
Intermediate mass, non-emission line B-type stars such as 
HR 2949 typically rotate even more slowly than this on the 
ZAMS, wi th their probability distributions peaking around 
Wo = 0.2 (iHuang et al.ll20lHl . 

Rather than try to estimate the maximum spindown 
time, we instead use the stellar age from Table [3] and 
the magnetospheric parameters from Table [6] to calculate 
IFo(iBp): 


Wo = 


(8) 


where fsp is the time before present in Myr. The results 
are shown in Fig. 1191 w here equatio n [5] is solved using rj 
as calculated using both IVink et al.l and IKrtickal mass-loss 
rates. _ 

Using Ivink et all ’s mass-loss rates, initially critical ro¬ 
tation is compatible with the extreme upper bound on HR 
2949’s age. If the age is closer to the nominal value of 11 
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Figure 19. The rotation parameter as a function of the time 
before present, tBP, for magnetosph eric parameters ca lculated 
with the mass-loss rates of IVink et alj (solid red line) an d iKrticd 
(dotted blue line). Vertical solid (green) lines indicate the stellar 
age, with the uncertainty indicated with dashed lines. Horizon- 
tal dash-dotted (blacl) lines Wq = 0.5, the theoretical maximum 
l|Lin et al.ll201lh . VF(tBp) is calculated such that fV(0) = W*. 
Dashed lin es indicate the uncertainty in W. In the case of the 
IVink et alJ mass-loss rates, 14^o = 1 can be acco modated within 
the uncertainty in the stellar age. Using iKrtickal mass-loss rates, 
the maximum value of Wq compatible with the stellar age is 0.55. 


Myr, however, equation [8] predic ts Wn c loser to the theo¬ 
retical upper limit of 0.5. Using iKrtickal ’s mass-loss rates, 
Wo < 0.5 at almost all plausible ages. 

The calculated spindown time, empirical and theoretical 
bounds on Wo, and the stellar age, are mutually consistent, 
and could be interpreted as indirect evidence for a magne¬ 
tosphere. However, the rotation parameter, W = 0.14lg J°, 
is also consi s tent w ith the probability peak observed by 
iHuang et al.l (l2010li . and hence with no spindown having 
occurred. Therefore, we conclude that while the stellar ro¬ 
tation is consistent with a magnetic spindown scenario, on 
its own this does not constitute evidence for the presence of 
a magnetosphere. 


9 CONCLUSIONS 

We have analyzed high-resolution spectropolarimetric 
observations of the optical double system HR 2948+2949. 
The principal results of this paper are that: 

1) HR 2949 is a member of the class of He-weak, magnetic 
Bp stars; 

2) Virtually all spectral lines, including H, are vari¬ 
able, and this variability furthermore distorts the magnetic 
diagnostics; 

3) There is no sign of a magnetospheric signature in 
the commonly used optical or UV diagnostic lines. 

Stellar parameters and rotational properties were de¬ 
rived from spectral line profile modeling, photometry, 
and evolutionary models. A probable misidentification 
of photometric measurements in the literature leads to 
different stellar parameters from those previously reported. 


with HR 2949 being hotter, more luminous, and larger than 
previously believed. 

Placement of HR 2948 on the Hertzsprung-Russell di¬ 
agram indicates it to be significantly older than HD 2949. 
We have considered different scenarios to explain this dis¬ 
crepancy, and conclude that the most likely are that either 
the binary system is not primordial, but formed via gravi¬ 
tational capture, or the stars are not physically associated. 

No magnetic field is detected in the single LSD measure¬ 
ment of HR 2948, with a longitudinal magnetic field error 
bar from the LSD profile of ~35 G. 

Zeeman signatures are detected in all of the 22 obser¬ 
vations of HR 2949. We find a rotational period of 1.90871 
d, and a dipolar model with i = 431^2° i = 50 ± 16°, 
and Bd = 2.4 ± 0.3 kG. This magnetic model is subject to 
revision, as the (Bz) curve itself is non-sinusoidal when mea¬ 
sured in both metallic lines and H lines. While a simple offset 
dipole model produces a better fit to Stokes V and (Bz) than 
a centred dipole, the fit is still unsatisfactory. Without per¬ 
forming Magnetic Doppler Imaging (MDI), departures from 
sinusoidal behaviour due to flux redistribution cannot be 
ruled out. The moderate inclination of the rotational axis, 
moderate usini, the large spectropolarimetric data-set, and 
the detectability of Zeeman splitting in numerous individual 
spectral lines, make HR 2949 a prime candidate for Doppler 
and Magnetic Doppler Imaging. MDI has been performed for 
only a small sample of Bp stars to date, and we still have 
a poor understanding of the prevalence of departures from 
dipolar magnetic fields in these stars, and whether these 
correlate in any way to e.g. the stellar age. 

Our abundance analysis indicates that HR 2949 is a 
He-weak Bp star, with significant overabundances of Fe- 
peak and rare earth elements, while HR 2948 is a chemically 
normal solar-abundance star. Precise abundances have been 
determined for only a very few Bp stars, making this anal¬ 
ysis a valuable contribution. Line-profile variability in He 
and metallic lines is almost certainly associated with photo- 
spheric chemical spots. Weak variability in the Balmer wings 
may be associated with variable displacement of H by other 
elements; variations in the strength of line blanketing; or 
a Lorentz force contribution to the photospheric pressure. 
Further modeling is required to distinguish between these 
possibilities. 

While there may be indirect evidence for the presence 
of a magnetosphere in the stellar rotation, as the maximum 
spindown age due to angular momentum loss via the mag¬ 
netized wind is in reasonable agreement with the age in¬ 
ferred from evolutionary tracks, there is no direct evidence 
for a magnetosphere, with H Balmer lines and UV resonance 
lines in absorption. Previously, this absence was something 
of a mystery, as the star’s magnetospheric parameters were 
quite similar to others with clear magnetospheric signatures. 
With our reanalysis of the stellar parameters, resulting in a 
substantially smaller Alfven radius, HR 2949’s character¬ 
istics are now within the range of other non-emission line 
magnetic B stars. Determining the conditions under which 
magnetically confined winds becomes visible requires stud¬ 
ies, such as the present one, that determine the stellar and 
magnetic properties of stars without emission, in addition to 
examination of stars with emission. Only then can the two 
groups be clearly separated in parameter space, an essential 
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condition for future clarification of the physics missing from 
current models. 
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Figure Al. As Fig.[7| for lines correlating with He. 







Figure A2. As Fig. [71 for lines exhibiting complex variations that neither correlate nor anti-correlate with He. 
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Table Bl. Longitudinal magnetic field measurements {Bz) and diagnostic null measurements {Nz) from LSD profiles. Measurements 
are given in order of acquistion; sets of measurements are ordered by decreasing number of lines included in the mask, given beneath the 
mask name. 



Metals 

Fe 


He 



S 



Si 



190 


50 


47 



40 



23 


Phase 

(Bz) (G) 

(Nz) (G) 

(Bz) (G) 

(Nz) (G) 

(Bz) (G) 

(Nz) (G) 

(Bz) (G) 

(Nz) (G) 

(Bz) (G) 

(Nz) (G) 

0.020 

446± 45 

-11± 45 

547± 89 

19± 89 

-190± 102 

-6± 102 

554± 

188 

-16± 

188 

610± 89 

78± 89 

0.076 

314± 30 

29± 30 

514± 56 

1± 56 

-219± 67 

-8± 67 

130± 

136 

15± 

136 

420± 55 

53± 55 

0.566 

-853± 29 

14± 29 

-769± 58 

5± 58 

-809± 47 

-23± 47 

-862± 

129 

-10± 

129 

-925± 58 

-57± 58 

1.000 

439± 26 

-4± 26 

672± 51 

-42± 51 

-219± 57 

9± 57 

307± 

110 

-9± 

110 

660± 53 

-23± 53 

0.082 

243± 28 

-31± 28 

375± 53 

-18± 53 

-275± 59 

5± 59 

198± 

130 

-115± 

130 

351± 51 

6± 51 

0.597 

-855± 31 

-5± 31 

-736± 64 

8± 64 

-760± 58 

10± 58 

-888± 

132 

7± 

132 

-919± 62 

45± 62 

0.094 

227± 28 

38± 28 

323± 55 

29± 55 

-283± 62 

21± 62 

430± 

125 

162± 

125 

421± 51 

-18± 51 

0.626 

-848± 28 

-25± 28 

-840± 55 

-93± 55 

-832± 45 

8± 45 

-891± 

123 

69± 

123 

-861± 57 

-19± 57 

0.178 

-95± 34 

49± 34 

82± 71 

14± 71 

-334± 64 

19± 64 

-74± 

179 

-133± 

179 

17± 62 

98± 62 

0.683 

-777± 39 

23± 39 

-764± 74 

103± 74 

-847± 58 

-29± 58 

-978± 

157 

-140± 

157 

-706± 81 

-24± 81 

0.156 

9± 33 

3± 33 

87± 66 

-76± 66 

-316± 74 

45± 74 

358± 

172 

-41± 

172 

66± 62 

-24± 62 

0.770 

-580± 41 

-38± 41 

-605± 76 

-81± 76 

-643± 61 

32± 61 

-665± 

164 

56± 

164 

-429± 92 

-12± 92 

0.288 

-489± 28 

42± 28 

-384± 66 

-33± 66 

-526± 47 

28± 47 

-428± 

150 

-69± 

150 

-516± 52 

-16± 52 

0.797 

-369± 32 

48± 32 

-226± 60 

17± 60 

-634± 53 

-34± 53 

-481± 

124 

88± 

124 

-227± 76 

5± 76 

0.348 

-630± 31 

31± 31 

-579± 75 

59± 75 

-605± 56 

17± 56 

-554± 

166 

-42± 

166 

-632± 59 

18± 59 

0.901 

179± 52 

40± 52 

407± 97 

145± 97 

-276± 86 

47± 86 

296± 

213 

27± 

213 

188± 113 

-196± 113 

0.489 

-712± 31 

-5± 31 

-645± 63 

66± 63 

-737± 52 

-21± 52 

-435± 

145 

93± 

145 

-810± 60 

11± 60 

0.370 

-636± 37 

-26± 37 

-682± 87 

80± 87 

-586± 64 

-78± 64 

-520± 

184 

-240± 

184 

-665± 70 

-64± 70 

0.510 

-866± 54 

24± 54 

-787± 112 

-58± 112 

-685± 87 

30± 87 

-568± 

248 

145± 

248 

-936± 101 

-31± 101 

0.411 

-739± 87 

91± 87 

-636± 186 

163± 186 

-654± 144 

10± 144 

-471± 

428 

-13± 

428 

-971± 168 

26± 168 

0.974 

400± 49 

56± 49 

764± 97 

91± 97 

-119± 101 

11± 101 

76± 

241 

143± 

241 

533± 102 

18± 102 

0.445 

-672± 22 

20± 22 

-667± 51 

70± 51 

-683± 38 

13± 38 

-709± 

131 

31± 

131 

-759± 46 

21± 46 


O 


G 


N 



Mg 



Ne 



21 


13 


10 



6 



3 


Phase 

(Bz) (G) 

(Nz) (G) 

(Bz) (G) 

(Nz) (G) 

(Bz) (G) 

(Nz) (G) 

(Bz) (G) 

(Nz) (G) 

(Bz) (G) 

(Nz) (G) 

0.020 

19± 158 

75± 158 

492± 135 

-90± 135 

1026± 293 

-125± 293 

47± 

217 

-21± 

217 

-338± 308 

232± 308 

0.076 

-76± 123 

31± 123 

344± 102 

-7± 102 

815± 201 

2± 201 

-59± 

134 

-110± 

134 

-196± 229 

-225± 229 

0.566 

-761± 136 

65± 136 

-1000± 85 

-59± 85 

-1188± 263 

239± 263 

-958± 

129 

-54± 

129 

-556± 155 

-170± 155 

1.000 

14± 97 

-117± 97 

491± 72 

123± 72 

807± 171 

169± 171 

89± 

133 

-128± 

133 

-460± 199 

-289± 199 

0.082 

-76± 116 

-96± 116 

211± 84 

-21± 84 

1136± 196 

121± 196 

-113± 

141 

-55± 

141 

-521± 210 

-131± 210 

0.597 

-631± 125 ■ 

■119± 125 

-1024± 101 

-41± 101 

-1398± 264 

84± 264 

-802± 

142 

-167± 

142 

-852± 188 

139± 188 

0.094 

-151± 107 

-22± 107 

139± 81 

-5± 81 

747± 189 

75± 189 

-211± 

146 

25± 

146 

-214± 218 

77± 218 

0.626 

-755± 133 

155± 133 

-1085± 77 

-11± 77 

-1275± 239 

-68± 239 

-812± 

138 

87± 

138 

-615± 175 

-116± 175 

0.178 

-382± 163 

216± 163 

-423± 105 

69± 105 

521± 289 

33± 289 

-537± 

158 

-19± 

158 

-596± 268 

215± 268 

0.683 

-690± 145 

-38± 145 

-1001± 110 

-31± 110 

-1382± 320 

281± 320 

-662± 

170 

-27± 

170 

-791± 194 

27± 194 

0.156 

-83± 135 

56± 135 

-91± 109 

17± 109 

487± 272 

-34± 272 

-501± 

160 

160± 

160 

-306± 245 

-401± 245 

0.770 

-564± 124 

-82± 124 

-689± 116 

28± 116 

-429± 318 

-80± 318 

-496± 

168 

-53± 

168 

-767± 190 

-97± 190 

0.288 

-622± 137 

-9± 137 

-658± 87 

-17± 87 

-427± 237 

13± 237 

-712± 

128 

147± 

128 

-525± 198 

129± 198 

0.797 

-489± 107 

-1± 107 

-407± 77 

75± 77 

-248± 207 

-35± 207 

-236± 

140 

41± 

140 

-558± 186 

18± 186 

0.348 

-657± 131 

14± 131 

-763± 97 

-11± 97 

-832± 290 

-203± 290 

-848± 

132 

86± 

132 

-619± 176 

-151± 176 

0.901 

-219± 150 

164± 150 

213± 146 

256± 146 

638± 340 

103± 340 

379± 

224 

111± 

224 

-503± 310 

19± 310 

0.489 

-708± 117 

20± 117 

-753± 89 

6± 89 

-756± 244 

-187± 244 

-865± 

134 

-152± 

134 

-709± 180 

96± 180 

0.370 

-739± 136 ■ 

■109± 136 

-880± 103 

-105± 103 

-484± 306 

143± 306 

-869± 

144 

18± 

144 

-464± 222 

184± 222 

0.510 

-949± 252 

-1± 252 

-989± 156 

-9± 156 

-707± 459 

10± 459 

-810± 

199 

-75± 

199 

-737± 335 

-83± 335 

0.411 

-756± 349 

345± 349 

-901± 267 

-300± 267 

-122± 585 

846± 585 

-456± 

326 

-215± 

326 

-539± 473 

-328± 473 

0.974 

-16± 127 

107± 127 

791± 162 

44± 162 

1014± 411 

734± 411 

-4± 

201 

123± 

201 

-149± 260 

220± 260 

0.445 

-677± 93 

-55± 93 

-861± 70 

-15± 70 

-895± 214 

-70± 214 

-756± 90 

79± 90 

-608± 99 

-10± 99 





22 M. Shultz et al. 


Table B2. Longitudinal magnetic field measurements (Bz) and diagnostic null measurements (Nz) from individual lines, arranged in 
order of increasing rest wavelength. 



Ca 11 

Si III 

Ti 11 


Fe II 


Hq 


-^0 

393.3663 

455.2622 

457.1968 

516.9000 

656.2797 

9 

1.170 

1.250 

1.572 


1.330 

1.000 

Phase 

0.020 

(B.) (G) 
101± 665 

{Nz) (G) 
129± 665 

{Bz) (G) 
759± 289 

{Nz) (G) 
113± 289 

{Bz) (G) 
168± 517 

{Nz) (G) 
386± 517 

{Bz) (G) 
726± 228 

{Nz) (G) 
-21± 228 

{Bz) (G) 
265± 46 

{Nz) (G) 
74± 46 

0.076 

69± 373 

-186± 373 

860± 154 

46± 154 

268± 449 

319± 449 

151± 185 

202± 185 

183± 29 

-51± 29 

0.566 

-864± 382 

-411± 382 

-793± 182 

-22± 182 

-771± 345 

-181± 345 

-949± 112 . 

■150± 112 

-733± 26 

-70± 26 

1.000 

-75± 310 

133± 310 

678± 160 

19± 160 

104± 339 

-333± 339 

448± 154 . 

■138± 154 

197± 31 

1± 31 

0.082 

163± 402 

552± 402 

440± 183 

105± 183 

1100± 459 - 

■1367± 459 

106± 137 ■ 

■190± 137 

116± 29 

-49± 29 

0.597 

-555± 354 

-77± 354 

-863± 183 

-198± 183 

-573± 463 

191± 463 

-742± 115 

82± 115 

-708± 41 

38± 41 

0.094 

-343± 341 

122± 341 

524± 191 

-260± 191 

350± 648 

239± 648 

-357± 175 ■ 

■172± 175 

0± 31 

7± 31 

0.626 

-1329± 370 

-230± 370 

-944± 201 

216± 201 

-58± 274 

-576± 274 

-823± 111 

-8± 111 

-661± 36 

-43± 36 

0.178 

-545± 456 

386± 456 

393± 155 

5± 155 

127± 499 

-6± 499 

60± 151 

164± 151 

-165± 41 

26± 41 

0.683 

-971± 476 

292± 476 

-687± 267 

66± 267 

-41± 434 

-272± 434 

-690± 158 

-38± 158 

-500± 41 

-24± 41 

0.156 

-161± 399 

-35± 399 

317± 169 

-298± 169 

264± 400 

-445± 400 

-91± 158 

33± 158 

-88± 36 

34± 36 

0.770 

-538± 530 

478± 530 

100± 302 

185± 302 

236± 441 

-340± 441 

-651± 164 

21± 164 

-258± 47 

-5± 47 

0.288 

-606± 355 

480± 355 

-567± 144 

49± 144 

-628± 361 

332± 361 

-422± 148 

-52± 148 

-513± 35 

17± 35 

0.797 

-847± 378 

-227± 378 

-29± 226 

-178± 226 

-570± 291 

225± 291 

-288± 107 

-46± 107 

-221± 31 

-2± 31 

0.348 

-696± 435 

-8± 435 

-549± 166 

107± 166 

-84± 455 

118± 455 

-241± 142 

148± 142 

-596± 35 

75± 35 

0.901 

247± 654 

-393± 654 

567± 268 

-200± 268 

159± 442 

-44± 442 

432± 247 

298± 247 

55± 52 

131± 52 

0.489 

-1116± 390 

-287± 390 

-753± 212 

-159± 212 

-540± 329 

400± 329 

-849± 105 ■ 

■346± 105 

-616± 32 

39± 32 

0.370 

-933± 439 

427± 439 

-927± 222 

-414± 222 

-1136± 555 

527± 555 

-732± 182 

152± 182 

-465± 43 

8± 43 

0.510 

-1054± 613 

332± 613 

-1012± 178 

11± 178 

-2105± 573 - 

■1449± 573 

-916± 187 ■ 

■763± 187 

-620± 32 

76± 32 

0.411 

-389± 979 

827± 979 

57± 199 

-25± 199 

-3686± 1551 1512± 1551 

-1336± 245 ■ 

-489± 245 

-582± 36 

-33± 36 

0.974 

284± 554 

-601± 554 

444± 218 

-158± 218 

94± 343 

350± 343 

1130± 124 . 

■318± 124 

34± 27 

21± 27 

0.445 

-773± 277 

-27± 277 

-500± 38 

-32± 38 

-603± 90 

-515± 90 

-737± 24 

98± 24 

-503± 6 

-22± 6 







